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ABSTRACT 

We revisit the paradigm of the dependency of jet power on black hole spin in accreting 
black hole systems. In a previous paper we showed that the luminosity of compact jets 
continuously launched due to accretion onto black holes in X-ray binaries (analogous 
to those that dominate the kinetic feedback from AGN) do not appear to correlate 
with reported black hole spin measurements. It is therefore unclear whether extraction 
of the black hole spin energy is the main driver powering compact jets from accreting 
black holes. Occasionally, black hole X-ray binaries produce discrete, transient (ballis- 
tic) jets for a brief time over accretion state changes. Here, we quantify the dependence 
of the power of these transient jets (adopting two methods to infer the jet power) on 
black hole spin, making use of all the available data in the current literature, which 
includes 12 BHs with both measured spin parameters and radio flares over the state 
transition. In several sources, regular, well-sampled radio monitoring has shown that 
the peak radio flux differs dramatically depending on the outburst (up to a factor 
of 1000) whereas the total power required to energise the flare may only differ by a 
factor < 4 between outbursts. The peak flux is determined by the total energy in the 
flare and the time over which it is radiated (which can vary considerably between out- 
bursts). Using a Bayesian fitting routine we rule out a statistically significant positive 
correlation between transient jet power measured using these methods, and current 
estimates of black hole spin. Even when selecting subsamples of the data that dis- 
regard some methods of black hole spin measurement or jet power measurement, no 
correlation is found in all cases. 

Key words: accretion, accretion discs, black hole physics, X-rays: binaries, ISM: jets 
and outflows 



1 INTRODUCTION 

Accretion onto black holes commonly results in relativis- 
tic, collimated outflows, or jets. In black hole X-ray bina- 
ries (BHXBs), these jets can brighten, fade or even switch 
off a ll together on timescales of days-months (see iFenderl 
l200rj . for a review). When jets are observed in BHXBs they 
usually manifest themselves as powerful, compact outflows 
commonly seen as a flat or slightly inverted radio spec- 
trum (a ~ - 0.5, wher e F v oc v a ) and are spatially 
resolved in some cases (e.g . iDhawan. Mirabel fc Rodrigued 
|2000| ; IStirling et "all l200ll . Miller- Jones et al. in prepa- 
ration). Both properties are in line with what is ex- 
pected for partially self-absorbed overlapping synchrotron 
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spectra |Blandford fc Konigll Il979l ). These compact jets 
of BHXBs are typically launched continuously for several 
months at least during a single outburst, when the X- 
ray spectrum is a hard powe r law (the 'ha rd state'; e.g. 
iMcClintock fc Remillardl |2006| ; iBellonil l20ld ) and are con- 
sidered to be analogous to persistent compact radio sources 
associated with active galactic nuclei (AGN). In BHXBs, 
most of the energy radiated by these jets is emitted at higher 
energies, for example in the infrared /optical regime ( e.g. 
ICorbel fc Fenderl l2002l : IChatv. Dubus fe Raichoorll201lf) or 
even higher - X-ray or possibly 7-ray e nergies in som e cases 



(e.g. iMarkoff Falcke fc Fended l200ll : iRussell et all l20ld : 



lLaurent et al.ll201ll ). The direct detection of a break in the 
jet spectrum in one BHXB, GX 339-4 at mid-infrared wave- 
lengths has allowed a measurement of the total compact 
jet radiative power to be made, > 6 x 10 35 erg s _1 (as- 
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suming a distance of 8 kpc) which i s > 0.6 per cent o f 
the total bolometric X-ray luminosity |Gandhi et al.ll201ll ). 
(and breaks have r ecently been constr a ined in several mor e 
BHXB jet spectra; iRahoui et ai1l2012l ; iRussell et al]|2013h . 
These jets are radiatively inefficient, and the total (ki- 
netic plus radiative) energy channelled into the jets can 
be approximately half of the avail a ble accretion energy 
or perhaps more jGallo et all 120051 ; IRussell et "ail l2007al ; 
iPakull, Soria fc Motchll20ld ). 

During accretion state transitions at high X-ray lumi- 
nosities (> 1 per cent of the Eddington luminosity; LBdd), 
a different type of jet is commonly seen in BHXBs. The 
compact jet fades and discrete, ballistic ejections appear 
as bright, optically thin radio flares, sometimes later re- 
solved into two expanding structures either side of the 
core BHXB (e.g. iMirabel fc Rodrigue j Il99i ; iTingav et ail 
1 19951 ). These discrete ejections are typically seen once, or 
at most a few times during a BHXB outburst cycle. Each 
ejection has a total radiat ive luminosity of ~ 10 ~ 



^Edd (see e.g. table 1 in iFender. Belloni fc GaUol |2004) . 
and may be associated with a fast, relat ivistic ejection (see 
also iMiller- Jones. Fender fc Nakarll2006l ). For comparison, a 
compact jet such as that of GX 339-4 mentioned above pro- 
duces this much ra diative energy in < 30 minutes during 
a bright hard state ([Gandhi et alj|2~01ll ). Since hard states 
normally last several months, the transient jets probably 
provide a negligible fraction of the total radiative power of 
jets from BHXBs. However, for both types of jet, the total 
kinetic power they contain is largely unknown, so compar- 
isons can only be made between their radiative luminosities. 
It is worth assessing whether the same physical process is 
responsible for launching both types of jet. 

The physical mechanism responsible for the launch- 
ing, acceleration and collimation of jets has been debated 
for decades, and remains one of the major paradigms 
in the field of accretion in strong gravitational fields. It 
is generally accepted that magnetohydrodynamic (MHD ) 
processes are at work (e.g. iMeier. Koide fc Uchidal l200lt) . 
but two groups of models are usually considered that 
can both produce the poloidal magnetic fields and 
differential rotation nece s sary to launch powerfu l jets 
|Blandford fc Znaiekl Il977l ; iBlandford fc Pavneill982l ). The 
Blandford-Znajek mechanism describes the extraction of 
the rotational energy of spinning BHs, in a region close 
to the BH where frame-dragging occurs (< 100 gravi- 
tational radii; R s — GM/c 2 ; e.g. IMeier. Koide fc Uchidal 



l200ll ; lMcKinnev, Tchekhovskov fc Blandfordll2012l ). The al- 
ternative to spin-powered jets is purely accretion-powered 
jets; the Blandford-Payne mechanism, where all the en- 
ergy and rotation originates in the disc and its mag- 
netic fields. The Blandford-Znajek model predicts that 
the total jet power scales with the second power of the 
dimensionless spin parameter, a* = cJ/GM 2 , where J 
and M are the black hole angular momentum and mass 
(see I Tchekhovskov. Naravan fc McKinnevI |2010| . for a sim- 
ilar scaling based on general-relativistic magnetohydrody- 
namics simulations) . 

The two leading methods currently adopted to infer BH 
spin in BHXBs are X-ray thermal continuum fitting, and X- 
ray reflection/Fe K a line fitting. The former applies to spec- 
tra of BHXBs in the thermal dominant (or h igh/soft) state 
l|McClintock fc Remillard]|2006l ; lBelloni|[201oh , while the lat- 



ter fits (low/)hard state spectra at Eddington rates higher 
than about 10 -3 . In essence, the thermal continuum fitting 
method provides a measure of the area subtended by the In- 
nermost Stable Circular Orbit (ISCO), while the reflection 
method returns a measure of the gravitational redshift that 
shapes the red wing of the broadened Fe K a line. Both meth- 
ods rely on the assumption that the accretion disc extends 
in to the ISCO, which i s in turn a fu nctio n of BH spin. The 
reader is referred to e.g. iMillerl (|2007h and lMcClintock et al.l 
(2011) for further details on the two methodologies. BHXBs, 
as opposed to AGN, are particularly suited for this kind of 
study. In outburst, Galactic BHXBs tend to be brighter than 
all AGN in the sky, so with more photons, spectral fitting 
can be achieved with a higher level of precision. In addition, 
the thermal continuum fitting method cannot be applied to 
AGN, since the peak of the disc emission falls in the UV 
band. 

Measures of BH jet powers also rely on different meth- 
ods, arguably all more indirect than those adopted for esti- 
mating the spin parameter, such as the relative normaliza- 
tion of the radio/infrared (IR) versus X-ray luminosity in 
the hard state for compact jets, and the total energy in the 
bright radio flares that are thought to be at the origin of the 
large scale, transient jets. 

In a previous work |Fender. Gallo fc" Russellll2O10l . here- 
after FGR10) we performed the first study to observation- 
ally test for empirical relations between BH spin (a*) and 
jet power. Radio-X-ray and IR-X-ray correlations exist in 
the hard state but some BHXBs are more radio-bright (or 
infrared-bright) than others, at a given X-ray luminosity. 
The luminosity of the jet at radio and IR frequencies rela- 
tive to the X-ray luminosity (the normalization of the cor- 
relation) was used to rank the BHXBs in order of jet power. 
No single relation was found between current reported BH 
spin measurements and relative luminosity of compact jets. 
This suggested that if the relative jet power estimates and 
BH spin estimates are accurate, compact jets may not be 
spin-powered. 

Since the acceleration and collimation regions of 
compact jets are usually consid ered to be at least 

tens to hundreds of R B (e .g. Junor. Biretta fc Livid 

1 19991; iMarkoff. Nowak fc Wilms! 120051; iGandhi et all l201lh 
bolko. Meier fc MarkoffeOloUPe'er fc Markofj|2012l '). these 
results may not be surprising because this is essentially out- 
side the region where BH spin energy can be tapped. In addi- 
tion, the jet is usually required to be anchored to a thick disc, 
which is thought to be present i n the har d X-ray state (e.g . 
lMeierlll999l ; IMarkoff et al.ll2005l ; see also lPolko et al.ll2010l ; 
however a thin disc that extends clos e to the ISCO may ex- 
ist ab ove 10~ Z/Edd in the hard state; iReis. Fabian fc Miller! 
|2010| ). Very recently, the actual base of the jet of the AGN in 
M87 has been spatially resolved using radi o interferometry, 
and i ts size was measured to be 5.5±0.47? g l|Doeleman et al.l 
l2012t ). However even here, the wide opening angle of the 
sub-relativistic jet and the existence of a counter jet indi- 
cates that this jet is anchored within (and powered by) the 
accretion disk and not lau nched due to magnetic field lines 
crossing the event horizon (jDoeleman et al.l 2012T ). 

The radio/I R normalization for a single source can 
change in time ([Russell et all l2007bl ; ICoriat et ail l201ll ; 
iGallo et all I2012T I which, if X-ray luminosity is not the 
quantity that is changing, invalidates a single relation- 
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ship between radio/IR luminosity normalization and jet 
power (or jet radiative efficiency). Some MHD simula- 
tions of disc-emitting jets have also successfully described 
prop erties of the outflows common l y see n in BHXBs 
(e.g. iFerreira et all l2006t iPolko et al l l2010h 



m 

In addition, 



iMigliari. Miller- Jones fc Russelll (|201lli found no significant 



correlation between the spin of neutron stars in X-ray bi- 
naries (which can be measured accurately for some but 
not all sources considered) and their relative jet powers 
inferred from radio or infrared jet luminosities (although 
less data were available for neutron star jets, and those 
data span a smaller range in luminosi ty than the s ample 
of BHXBs). Finally, in a related study. iBanibl HoH) con- 
siders the Johannsen-Psaltis metric around a BH, a specu- 
lative, alternative space-time metric to the Kerr metric. In 
this scenario, current data of BHXBs suggest there may be 
a correlation between BH spin as defined by this metric, and 
compact jet power. 

The only tentative relation found by FGR10 was be- 
tween transient jet ejection energy and BH spin mea- 
sured via reflection fits (but not continuum fits), from 
four data points (fig . 6 in FGR10). More recently, 



iNaravan fc McClintocM (2012; hereafter NM12) claim a re- 
lation between transient jet peak radio luminosity and BH 
spin measured from thermal continuum fitting, using a sub- 
sample of the data in the literature. NM12 take the peak 
luminosity of the radio flare associated with the hard-to- 
soft state transition as their estimate of the jet power (as 
opposed to an estimate of the total energy of the ejection). 

Here, we test for a relation between transient jet power 
and reported BH spin using all the available data (de- 
scribed in Section 2), adopting two differing methods to in- 
fer the jet power and two methods to estimate BH spin. We 
make use of a Bayesian fitting routine in order to quan- 
tify the significance of any relation (the methodology is 
explained in Section 3). In all cases, no significant cor- 
relation is found. If the methods used to infer jet power 
and BH spin are correct, we can rule out a single posi- 
tive correlation between transient jet power and BH spin 
with high confidence. A discussion and brief summary is 
provid ed in Section 4. We also comment on a very recent 
paper, ISteiner. McClintock fc Naravanl (|2013h in Section 4. 
This new work follows on from that of NM12, adding one 
more source to their sample. This new source was already 
included in our samples here, and we demonstrate that this 
paper does not affect our results or conclusions. 



2 THE SAMPLE 

In Table 1 we list all BHXBs that we are aware of with both 
BH spin (measured either by the continuum fitting method 
or the reflection fitting method) and observed transient jet 
flares. BH spin a nd radio data are t aken from table 1 of 
FGR10, table 1 of lFender et alj l|2004l ) and recent estimates 
noted in NM12, plus some new BH spin estimates published 
since FGR10 but not used in NM12 and radio flare data from 
the same BHXBs (see Table 1 for all original references for 
the BH spin and radio data, as well as distances and BH 
mass estimates). Where more than one spin estimate exists 
using the same method for the same BH, we use the most 
recent spin measurement, inferred using data from modern 



X-ray telescopes, since these are likely to invoke the most 
sophisticated fitting techniques and use the highest quality 
data. 



2.1 BH spin parameters 

Throughout this paper we consider both leading methods 
of measuring BH spin, and do not favour one method over 
the other. It is not in the scope of this paper to argue that 
one method provides more accurate spin estimates than the 
other; a scientifically fair study should consider all possible 
published measures. In the long term, the hope is that the 
two methods will converge and give consistent results. This 
appears to be currently the case for th ree BHXBs (GRS 
1915+105, Cyg X-l and XTE J1550-564; [McClintock et all 
| 2006l; iBlum et al.|[2009l ; ISteiner et al.ll201ll ; iGou et al.ll201lT 
iFabian et al.ll2012l ). 

For some BHXBs, multiple spin measurements us- 



ing the same method have been broadly consistent with 
each other, whereas for other BHXBs they have not. 
Two examples are the various spin estimates based on 
the disc continuum fitting meth od for GRO J 1655-40; 

1997h- 



0.93 



a* 



0.78 ± 0.10 



jGierlinski. Maciolek-Nie'd zawiccki & Ebisawa 2 00l|). a* < 
0.7 |Davis fc Hubenvll2006l ). a, = 0.7 ± 0.1 jshafee et all 
l2006t ), and t he reflection fittin g method for Cy g X-l; a« = 
0.05 ±0.01 (|Miller et al.ll2009l). a. = 0.88±£ ?I (|Duro et all 



1201 J ), a, = 0.97+Q 02 4 (|Fabian et all \20m . In almost all 



cases the most recent measurement is likely to be the most 
accurate one, but we caution that this may not always be the 
case, and it depends on the many subtleties of the techniques 
and assumptions adopted in each individual case, which we 
cannot go into here. Generally speaking, for most BHXBs 
the error bars of most, if not all measurements overlap when 
one spin measurement method is considered, but in a few 
cases there are large differences (like that of Cyg X-l men- 
tioned above; see also the discussion in FGR10). 



2.2 Jet power from radio flares 

Since the peak flux of the bright radio flare of a BHXB 
can differ dramatically for different outbursts of the 
same source, we tabulate the brightest radio flare ob- 
served during state transition for each outburst of each 
source. In addition to the sources considered here, there 
are several other notable cases of XBs with a large 
range in peak radio flare fluxes (e.g. Cyg X- 3 , Cir 



X-l, SS 433; iPreston et all [l983l; Ispencer et al.l 19861 ; 



iBonsignori-Facondi etal.1 1 19861; iMioduszewskjetalJ 200 ll; 



Tudose et al. 2008; Koljoncn et al. 
201 it iBlundell Schmidtobreick 



Corbel et al.ll2012h . 



120101; IWilliams et all 
fc Trushkinl l201ll ; 



For one source, 4U 1543-47, a radio flare seen during 
the state transition during its 2002 outburst was one of the 
few in the literature wher e the rise of the fl are was reported 
and well sampled in time l|Park et alj2004h . rising at 0.8-1.0 
GHz from 6 ± 2 mjy to 16 ± 2 mjy in 6.0 h, then peaking 
at 21.9 ± 0.6 mjy just 2.0 h later. Less than a day (19.0 h) 
after this brightest detection, the source had faded to < 2.7 
mjy at 0.8 GHz. The light curve morphology is typical of 
a flare event often seen over the state transition (see e.g. 
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Table 1. All BHXBs with both published radio jet flare data and BH spin estimates. 



Source 


d 




BH spir 


a, — 


^5GHz 


At 


Year 






References 






(kpc) 


(M e ) 


disc 


reflection 


(mjy) 


(h) 




d 


M BH a, 


•S'5GHz 


A0620-00 


1.06 ± 0.12 


6.61 ± 0.25 


0.12 ± 0.19 




203 




1975 


1 


1 


2 


3,4 


GRS 1124-68 


5.1 


6.0 + 1.5 


-0.04° 




171 




1991 


5 


6 


7 


8 


4U 1543-47 


7.5 ± 1.0 


9.4 + 1.0 


0.8 + 0.1 


0.3 + 0.1 


17.5 


8 


2002 


9 


9 


10,11 


12 


XTE J1550-564 


4.38 + 0.50 


9.10 ± 0.61 


0.34 + 0.24 


<^+ - 15 

U.OO_ Q 22 


265 


12 


1998 


13 


13 


14 


15,4 


XTE J1652-453 


8 6 


10 b 




~ 0.5 


> 0.95 




2009 


_6 


_b 


16 


17 


bKU J 1655-40 


3.2 ± 0.5 


c on I n T7 

6.30 ± 0.27 


0.7 ±0.1 


n < \ o i n r\ i 

0.98 ± 0.01 


2 420 
3.6 


12 
71 


1994 
2005 


18 


19 


10,11 


18,4 
20 


GX 339-4 


8 b 


12.3 ± 1.4 


< 0.9 


0.94 + 0.02 


55 


5.5 


2002 


21 f ' 


22 


23,11 


21 


H1743-322 


8.5 ± 0.8 


13.3 ± 3.2 


0.2 + 0.3 




96.1 
24.4 


<48 
<24 


2003 
2009 


25 


22 


25 


26 
27 


XTE J1752-223 


3.5 ± 0.4 


9.5 + 0.9 




52+ 016 


20 




2010 


28 


28 


29 


30 


GRS 1915+105 


11.0 + 1.0 


14.0 + 4.4 


0.975 ± 0.025 


0.98 + 0.01 


912 
356 
150 
356 


<119 
12 

<24 


1994 
1997 
2006 
2010 


31 


32 


33,34 


35,4 
31,4 
36 
37 


Cyg X-l 


Last?:!; 


14.8 ± 1.0 


0.96 + 0.04 


97+0-014 
u,a -0.02 


186 


>0.3 


2005 


38 


39 


40,41 


42 


GS 2000+25 


2.7 + 0.7 


7.2 + 1.7 


0.03° 




> 4.9 


<525 


1988 


43 


43 


7 


41 



a For GRS 1124-68 and GS 2000+25, BH spin estimates were reported with no errors quoted; here we assum e the same errors a s 
A0620— 00 (which has a similarly low spin). 6 The distance and BH mass of some BHXBs are uncertain (e.g. iHvnes et al]l2004l ; 
iHiemstra et alj|201ll ); here we assume 8 kpc (since these sources are towards the galactic centre) and 10 Mq, respectively. The peak 
radio flux density and rise times are tabulated for each radio flare over the hard-to-soft state t ransition reported for eac h source. In the 
sixth column, the year is t hat of the radio flare. Reference s : (1) = ICantrell et alj ll20ld); (2) =lGou et alj l l2010h; (3) =lKuuIkers et al.l 
lll999l); (4) = NM12; (5) = lGelino, Harrison fc McNamaral j200lT): ( 6) = Remillard fc McCIintockl ll200d) ; (7) = lzhang et alj 1119971 1; (8) 
= lBall et alj j 19951) ; (9) = lQrosz et alj i2002l); ( 10) = IShafee et alj j2009): ( 1 1) =lMiller et alj J2009h; (12) =|Park et alj ( 20041); (13) = 
lOrosz et alJll2011bl): (14) =ISteiner et alj I^OllI ); (15) =lHannikainen et £d1 ||2009| ): (16) =lffiemstra et al.l ll201ll); (17) =ICa,Ivelo et alj 
||2009|); (18) =lHiellming fc Rupenl lll995l); (19) =lGreene Bailvn fc Orosj j200ll); (20) = lRupen Dhawan fc Mioduszewskil d2005T) ; (21) 
= lHvnes et all <2004l): (22) = IShaposhnikov fc Titarchukl <2009l) ; (23) = iKolehmainen fc Pond foOld); (2 4) = lGallo et alj ||2004|): (25) 
= ISteiner, McCIintock fc Reidl j2012T ); (26) =lMcCIintock et all J2009I); (27) = [ Miller- Jones et al.N2012h; (28) = IShaposhnikov et ail 
i2O10D ; (2 9) =lReis et alj <201ll); ( 30) = lYang et alj feOld); ( 31) = IFender et alj lll999l) ; (32) = lGreiner, Cubv fc McCaughreanl d200lT) ; 



(33) =lMcClintock et alj 1120061) ; (34) _ .. 

iTrushkin fc Nizhelskiil |201C|); (38) =[Reid et alj feOl 
l20lj; (42) 



Rodrig;uez et alj jl995l); (36) =lRushton et al.1 ll20ld); (37) 
Orosz et alj ll2011al) ; (40) = lGou et alj ll201ll); (41) = iFabian et alj 



Blum et alj d2009fl ; (35, 
. i et alj l l201ll); (39) =|Prosz et alj ll2011al) ; 
Fender et akl I2OO6I) ; (43) = lBarret, McCIintock fc Grindlavl jl99d) ; (44) = lHiellming et al.1 l| 19881) 



iGallo et al.|[20o3 ) . Radio flares typically rise faster than they 
fade, so it is unlikely that the peak of this radio flare was 
much brighter than the observed peak of 22 mjy (NM12 
take this flux density as a lower limit). In addition, the rise 
times of all radio flares are consistent with being aroun d 
5 - 12 h (Table 1; see also table 1 of iFender et~al1 l2004 l. 
As an extremely conservative estimate, the peak could have 
been twice as bright as that observed, so we estimate a peak 
radio flux density of 22-44 mjy at 0.8-1.0 GHz. We use the 
method of NM12 of assuming a = —0.4 for the spectral 
index of the optically thin radio flare in order to estimate 
the 5 GHz flux, which was 11.6-23.4 mjy. 

GRO J1655-40 had a very bright radio flare of 2 
Jy in 1994 which was associated with relat ivistic ejecta 
ijHiellming fc Rupenl [l995l ; iTingav et al. I ll995l) . whereas in 
2005, an optically thin radio flare was detected at 6.6 ± 0.4 
mjy, 3.6 ± 0.1 mjy and 3.8 ± 0.3 mjy at 1.43 GHz, 
4.86 GHz and 8.46 GHz, r espectively, which then decayed 
over the next few days |Rupen et al.1 120051 ). The radio 
light curve is well sampl e d, with daily coverage (see also 
IShaposhnikov et al.1 120071 ; IFender. Homan fc Bellonl 120091 . 
and references therein) and the rise and decay of the flare 
were observed. This was the brightest flare seen from GRO 
J1655-40 during the hard to soft state transition of this 
outburst, and occurred within days of the source entering 



the soft state (|Rupen et al 1120051: IShaposhnikov et al.ll2007l ). 
which is fairly typical ( Fender et alj|2009l ). From these two 
outbursts with well sampled radio monitoring, it is clear 
that the peak luminosity of the flare associated with the 
state transition can vary by three orders of magnitude. The 
fainter flare in 2005 had a longer rise time, and the resulting 
total energy required to produce the 1994 and 2005 radio 
flares may only differ by only a factor ~ 4 (see Section 3). 

Four bright 5 GHz radio flares of GRS 1915+105 are 
known to have peaked at flux densities between 150 and 912 
mjy (Table 1). NM12 already noted that the bright radio 
flare of GRS 1915+105 seen in 1997 peaked at a luminosity 
~ 3 times fainter than the bright radio flare in 1994, which 
invalidates a single relation between BH spin and peak radio 
flare luminosity because the BH spin cannot change in three 
years (NM12 only included the brightest of these two flares 
of GRS 1915+105 in their following analys is of the BH spin- 
jet power relation). iRushton et al. (|2010l ) found from high 
resolution radio imaging of GRS 1915+105 that even fainter 
radio flares are resolved into relativistic discrete ejecta. We 
therefore include all four bright radio flares in Table 1 and 
the following analysis. These bright radio flares are not to 
be confused with the fainter, quasi-regular 'oscill ation' flares 
which peak at ~ 50 mjy (e.g. IFender et al]l2004l ). and which 
are not included here (but which could be due to the same 



© 2012 RAS, MNRAS 000,[THTT] 



On the powering mechanism of transient jets 5 



(a): BH spin from disc, Pj et from radio flare total energy: 
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(b): BH spin from reflection, P| et from radio flare total energy: 
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(c): BH spin from disc, P; et from radio flare peak luminosity: 
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(d): BH spin from reflection, Pj et from radio flare peak luminosity: 
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Figure 1. BH spin versus relative jet power. In the left panels BH spin measured via the disc continuum method are plotted; in the 
right panel BH spin measured via the reflection method are plotted. The jet powers es timated from total e nergy (FGR10) and peak radio 
luminosity (NM12) and shown in the upper and lower panels, respectively. NM12 and lSteiner et ah! (120131 ) used a subsample of the data 
shown in panel (c). 



physical process). HI 743-322 also had peak radio flare flux 
densities in 2003 and 2009 that differed by a fac tor of ~ 4 
(|McClintock et al.|[2009l ; iMiller- Jones et al.ll2012t ). 

We also include Cygnus X-l in our analysis. NM12 dis- 
card this source because, as they argue, the mass transfer 
is via a stellar wind from its high-mass companion star. Al- 
though Cyg X-l is a wind-fed system, an accretion disc has 
been detected around the BH, its size has been estimated 
dCoriat. Fender fc Dubusll2012(h and the disc has even been 
used to infer the BH spin fe.g. lMiller et al.ll2009l ;l Gou et al.l 
Eon]). On large scales, the mode of accretion for wind-fed 
systems differs to that of Roche lobe overflow systems - the 
stellar wind ensures a relatively steady, high mass accretion 
rate onto the outer accretion disc. The influence of BH spin 



on jet production is only important on distances of <C 100i? g 
from the BH in which general relativisti c frame dragging 
becom es important (e.g. lMeier et al.ll200ll ; [McKinnev et al.l 
2012). Obviously the mode of accretion towards the BH on 
these much smaller size scales is the same for wind-fed sys- 
tems as (transient) Roche lobe accreting systems (i.e., mat- 
ter is accreted from the inner disc into the inner regions that 
are affected by frame dragging), and is not affected by the 
mode of accretion on large scales (accretion disc or wind-fed 
accretion). 

In addition, it is well known that Cyg X-l performs 
state transitions that are associated with radio flares (e.g . 
iFender et aTl 120061 : IWilms et aLll2007l ; iRushton et~aH 120121 ) 
like other BHXBs. The brightest radio detection of Cyg X-l 
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ever seen over years of detailed, long-term monitoring with 
the Ryle and Green Bank te lescopes (e.g. iBrocksopp et al] 
1 19991 ; lLachowicz et al] 120061 ) was during the hard to soft 
state transitio n in 2004, where th e radio luminosity peaked 
at 120 mJy (jFender et al] 120061 , in addition, a likely ex- 
tended, discrete ejection has been seen). In fact Cyg X- 
1 spends 10 per cent of its time in the soft state (e.g. 



iGallo et al. 20051) in which the radio emission is suppressed 
(e.g. iTigelaar et al][20ol ; iRushton et al.1120121 ). as expected 
after crossing t he 'jet line' when th e ballistic ejection is usu- 
ally launched (jFender et al] 120041 ). Cyg X-l should there- 
fore certainly be included; here we use the aforementioned 
brightest radio flare of 120 mJy at 15 GHz. Again adopting 
a = -0.4 (NM12) we estimate a 5 GHz flux density of 186 
mJy. 



3 ANALYSIS AND RESULTS 

The two methods we adopt to estimate the relative jet power 
of the discrete ejections are explained in FGR10 and NM12, 
respectively. In FGR10 we adopted the classical calculation 
of the minimum total energy required to energise a syn- 
chrotron fla re, P oc (cAt) 9 ^ 7 L 4 ^ 7 where L is the peak radio 
luminosity |Burbidgelll959l ; lLongai n ll994l ). This calculation 
is frequently used in the study of synchrotron radiation from 
numerous sources at wavelengths from radio to X-ray, and 
is a reliable measure the minimum total luminosity of an 
individual flare. We calculate the jet power of each radio 
flare listed in Table 1 as a fraction of the Eddington lumi- 
nosity (using the known values of distance and BH mass). 
The second method of measuring jet power is simply using 
the luminosity at the peak of each radio flare (NM12), which 
assumes the peak luminosity is a linear indicator of total jet 
kinetic pow er (contrary to the generally accepted theoretical 
work s; e.g. iBlandford fc Konigll Il979l: iFalcke &i Biermannl 



1 19951 ; iMarkoff et al] 120031 ; iHeinz fc Sunvaevl 120031 ). This 
method contains fewer uncertainties, but does not take into 
account the different temporal profiles and hence different 
total energies of each radio flare. Note that for both meth- 
ods, the jet power as a fraction of the X-ray luminosity is 
not considered, because quasi-simultaneous X-ray luminosi- 
ties were not available for all radio jet flares. 

In Fig. 1 we plot BH spin versus relative jet power for 
all data listed in Table 1. We assume an error of a factor of 
two for the jet power estimates. For those sources for which 
different measurements/limits of the jet power have been re- 
ported over different epochs, the highest values correspond 
to filled black markers, while grey markers correspond to 
lower values/limits. It is clear that there is a large scatter 
in the values of jet power for a single source, because tran- 
sient jets had different peak fluxes and total energies during 
different outbursts. 

BH spins of a, 0.7 measured from disc continuum 
fitting have a range of jet flare powers spanning ^ 2.5 orders 
of magnitude (left panels of Fig. 1). With the inclusion of all 
the available data, the claimed Pj e t oc relation of NM12 is 
no longer as obvious (Fig. lc). Particular sources to note are 
4U 1543-47, Cyg X-l and GRO J 1655-40 (2005 outburst), 
which all had radio flares orders of magnitude fainter than 
expected from the relation, whereas the radio flare of GRS 
1124-68 was two orders of magnitude too bright for its low 





Figure 2. Linear regression analysis results for spin parameter as 
estimated from thermal continuum vs. jet power from the peak lu- 
minosity of the radio flare (after NM12). The right column shows 
the posterior distribution of the fitted slope, with negative values 
highlighted in black. The left column shows a random selection 
of 150 out of the 10,000 fitting iterations that are represented to 
the right, over-plotted to the data with grey dotted lines. The top 
panels include the full data set discussed in this paper; the mid- 
dle panels treat a simulated distribution of data points with the 
same values for the spin parameter as in the top panels, but jet 
powers randomly drawn from a positive correlation of the form 
Pjet oc a\ . The bottom panels repeat the same exercise as in the 
top panels, but for an arbitrary sub-sample of data points. The 
labels in the right column denote the inferred range of values for 
the slope of the correlation, while the percentages refer to positive 
slopes. Percentages lower than ~ 95 correspond to no statistically 
significant positive correlation (since more than ~ 5 per cent of 
the draws result in a negative slope). 



BH spin (although we caution the reader that this is a BH 
spin measurement inferred from old data, which as yet has 
not been superseded). 

Spins of a, ^ 0.9 measured via reflection have jet pow- 
ers ranging ^ 1.5 and ~ 3 orders of magnitude for jet powers 
estimated from their total energies and peak luminosities, 
respectively. There is no clear single relation between jet 
power and BH spin in any of the four panels in Fig. 1. 



3.1 Statistical analysis 

In order to quantify our visual assessment, we investigate 
the presence of a positive correlati on between jet power 
and BH spin using a Bayesian code (|Kellvll2007h . We carry 
out a linear regression analysis on a relation of the form 
log (Pjet) — a + /31og(a*), with intrinsic random scatter in- 
cluded, and for the four possible permutations between the 
two different methods for estimating the transient jet power 
and and spin parameter, respectively. Three Gaussians (for 
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the slope, intercept and intrinsic scatter) are used in the in- 
dependent variable mixture modelling, and a minimum of 
10,000 iterations are performed. The most likely parameter 
values are estimated as the median of 10,000 draws from the 
posterior distribution. For those sources for which different 
values of the radio flare luminosity have been reported over 
different epochs, we adopt the highest value as the tracer for 
jet power (corresponding to the filled black marks in Fig. 1). 

The top panels of Fig. [2] illustrate the results of the 
linear regression analysis for the case where a correlation 
has been claimed, i.e. for the jet power as inferred from the 
peak luminosity of the radio flare versus the spin parame- 
ter as estimated from the thermal continuum fitting method 
(NM12), here denoted as ad. The histogram on the top right 
shows the distribution of the fitted slopes, with the draws 
returning a negative slope highlighted in black. For this spe- 
cific case, only 60 per cent of the draws have slopes consistent 
with a positive correlation between the spin parameter and 
and jet power being present, ruling out a positive correlation 
with high confidence (> 5 tr). This is perhaps most apparent 
in the top left panel, where a random sub-sample of 150 fits 
(out of 10,000) is over-plotted to the data. 

For comparison, the middle panels correspond to a sim- 
ilar analysis as the one discussed above, only for a simulated 
data set where the spin parameter values (along with errors) 
are identical to those in the top panels, while the jet power 
values are randomly drawn from a Pj e t oc o» distribution 
with an intrinsic scatter of 0.5 dex (additionally, each jet 
power value is assigned an error of 0.3 dex). In this case, 
only 0.4 per cent of the draws return negative slopes, which 
indicates a statistically significant (> 3 a) positive correla- 
tion. By comparison with the top panels, this exercise serves 
to illustrate that the lack of a statistically significant corre- 
lation in the full data set is not merely a result of the small 
number statistics, and the code could indeed recover a sta- 
tistically significant correlation with eight data points only. 

Lastly, the bottom panels illustrate how the results 
of the regression analysis change for five out of the eight 
data points considered in the top panels (specifically with 
GRS1124-68, Cyg X-l and 4U 1543-47 excluded from the 
analysi s; as in NM12). Th e sample is actually the same as 
that in ISteiner et af1l|2013h ; this is discussed further in Sec- 
tion 4. While the slope distribution is formally consistent 
with a value of two, the significant tail of negative draws 
(11 per cent) corresponds to a marginally significant corre- 
lation only (< 2 <t), although it is worth stressing that the 
statistical significance of the correlation for this data set is 
very sensitive to the error bars (the Bayesian code we are 
adopting can only treat with symmetric error bars. 

The same analysis as described above was carried out 
to assess the presence of a relation between the spin pa- 
rameter a and the jet power p as measured from the differ- 
ent methods, namely spin from reflection (o r ) and jet power 
from total energy (denoted as Pfgr, vs. Pnm). Plots analo- 
gous to those discussed in Fig. [2] are shown in Fig. [3] for the 
following combinations of parameters: jet power from peak 
radio flare luminosity vs. spin parameter from thermal con- 
tinuum (top left); jet power from peak radio flare luminosity 
vs. spin parameter from reflection (bottom right); jet power 
from flare total energy vs. spin parameter from thermal con- 
tinuum (top left); jet power from flare total energy vs. spin 
parameter from reflection (top right). As apparent from the 



prominent negative tails of the fitted slope distributions, no 
positive, statistically significant correlation is found in all 
cases. A similar conclusion is reached by fitting either p as 
a function of a (solid lines), or a as a function of p (dashed 
lines) . 



4 DISCUSSION AND SUMMARY 

We have shown that no single statistically significant re- 
lation exists between inferred transient jet power and BH 
spin parameter, with the current estimates of these param- 
eters. Furthermore, the very fact that one BHXB can pro- 
duce radio ejections with very different total energies and 
peak luminosities invalidates a single relation between tran- 
sient jet radiative energy and BH spin. This latter discovery 
implies that many jet flares may be less luminous than the 
maximum value that would be produced by a purely spin- 
powered jet. For a typical source, there have been a lim- 
ited number of outbursts observed to date, so it could be 
that truly spin-powered jets do reach a constant, maximum 
jet power, but few flares reach this critical maximum. This 
could lead to underestimated jet luminosity estimates. We 
note that one source, GRS 1124-68 appears to have a high 
jet power but a low BH spin, which is not consistent with 
the above picture. Regular radio monitoring of populations 
of BHXBs over state transitions is required to measure the 
distribution of jet powers of discrete ejections. 

It is worth noting that in our study we do not include 
the effects of relativistic beaming of the jet emission. This is 
because both the orbital inclination angles (which are usu- 
ally assumed to be perpendicular to the jet axis) and the 
jet Lorentz factors are uncertain for many sources (NM12 
assumed a Lorentz factor of 7 = 2 for all sources). 

Very recently, a new paper has appeared from the same 
group as NM12, which again claims a correlation between 
BH spin (measured from continuum fitting) and jet power 
(measured from pe ak luminosities) in transient ejections 
l|Steiner et al.ll2013h . In this work, the same subsample of 
BHXB spin measurements that appeared in NM12 was used, 
with an addition of one more spin measurement, that of 
H1743-322 from ISteiner et al.1 (|2012l ). This spin measure- 
ment is also included in our analysis. However, it is worth 
noting a few discrepancies between our analysis and theirs. 

For this source we to ok jet flares f r om tw o outbursts, 
2003 and 2009, whereas in lSteiner et all i|2013l ') just the for- 
mer jet flare is used. Furthermore, we adopted the radio 
flare peak luminosity associated with the hard to soft tran- 
siti on in Fig, lc (the sam e method as was adopted in NM12) 
but ISteiner et al.1 l|2013h use a fainter flare, which peaked at 
35 mjy at 5 GHz instead of the brighter, 96 mjy flare we 
used. Their arguement is that this fainter flare must have 
been the one associated with the launch of the jet ejection 
because the date on which this ejection occurred is known 
from resolved radio observations. However, if the BH spin 
powers the 35 mjy flare but not the 96 mjy flare, then 
what powers the brighter flare? The 96 mjy flare had an 
optic ally thin spectrum and did not occur during the hard 
state (|McClintock et al.ll2009h . so we can rule out a compact 
jet origin to this brighter flare. If one of those two flares was 
associated with an energetic boost from the spin of the BH, 
one would expect it to be the brighter one. 
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Figure 3. Same as Fig. [2] but for different combinations of methods for estimating the spin parameter a and jet power p. a T and a<j 
denote spin parameters from reflection and disc fitting, respectively. Pfgr and Pimm refer to jet power from total energy and peak flare 
luminosity, respectively. The solid/dashed lines are the fitted p(a), a(p) slopes, compared to a p oc a 2 slope (dotted line). Dotted lines in 
the histograms are for slopes of 2 (p vs. a fit) and 0.5 (a vs. p). The numbers in the two histograms next to each fit are the percentage 
of draws that have slopes consistent with some positive correlation between a and p being present (i.e. the unfilled fraction). The data 
points do not show evidence for a significant correlation, for any of the four a versus p combinations. 



Furthermore, since we showed that radio ejections ob- 
servationally possess very different peak luminosities, the 
peak flux cannot be used as a proxy for the total power 
of a discrete ejection. ISteiner et al.l (|2013h do not take into 
account the time over which a flare is radiated (which empir- 
ically can vary by more than one order of magnitude; see Ta- 
ble 1), but instead adopt [vL v ) T ^ oc P\ et (see equation B10 
in the Appendix of lSteiner et al.f201St ) , contrary to standard 
classical synchrotron equation, Pj c t,min oc (cAt) 9 ^ 7 L 4 ^ 7 . 

Finally, a further discrepancy lies in the assumed rela- 
tion b etween BH sp i n and jet power. In fitting the corre- 
lation, ISteiner et atl l|2015t) assume Pj c t oc (a„/(2(l + (1 — 
a *) 1//2 ))) 2 ! as opposed to that adopted in NM12, Pj ct oc al, 
assume the same jet Lorentz factor for all sources, and fix 
the power law index to 2.0 (as indicated in footnote 6 in 
their paper). For comparison, the Bayesian, linear regression 
analysis we present here leaves the slope as a free parameter, 
testing the probability that it is positive at all. As we dis- 
cussed in detail in Section 3, a full Bayesian treatment, with 
error bars on both jet power and spin parameters included, 
shows that a positive correlation of the form Pj ct oc al is 
at best marginally significant (< 2 a), even when the new 
H1743-322 data point is included. 

If any of the above methods of measuring BH spin or 
jet power are correct, then we can conclude that there is 
no single relation between BH spin and transient jet power. 
Are there any reasonable explanations that could explain 
the lack of a correlation? Many works predict (analytically 



or via simulations) a relation between the spin of a black 
hole and the power of the r esulting accretion-fed jet (e.g . 



MacDonald fc Thorndll982i; iLivio Ogilyie fe Prin glc 1999; 
McKinnevi 120051; |Pe Villiers et al.l 120051; lHawlev fc Krolikl 
20061; iTchekhovskov. Naravan fc McKinnevi |2010| ; 
McKinnev, Tchekhovskov fc Blandfordl l2012i ). Positive 



relations are generally predicted but not necessarily as 
simple as the Pj e t oc al relation originally derived by 
iBlandford fc Znaiekl (|l977T I. It has recently been shown 
that rapidly retrograde spinning BHs (where the BH 
spin is in the opposite direction to the accreting matter 
spiralling inwards) can produce jet s just as powe r ful as 



Garofald [20091 : iMeierl 
20121 ). 



Interestingly, 



rapid ly prograde spinning BHs (e.g . 
l201ll : ITchekhovskov fc McKinnev 
the only source with a (slightly) negative spin in Table 
1, GRS 1124-68, had a relatively powerful jet (Fig. 1; 
lower left panel). However, since several BHs with very 
similar spin estimates have produced flares of energies and 
peak luminosities spanning several orders of magnitude 
(including different jet powers for different flares from the 
same source) we can rule out a single relation between jet 
power of discrete ejections and current estimates of BH 
spin. This is quantified by our Bayesian analysis (Figs. 2 
and 3). 

The spin parameters as estimated from the two methods 
do not agree with each other for a number of sources, and 
neither is to be preferred. In addition, for compact, steady 
jets that exist in the hard state, radio or infrared luminosity 



© 2012 RAS, MNRAS 000,[THTT] 



On the powering mechanism of transient jets 9 



normalization may turn out to be an unreliable tracer for jet 
power. Dual tracks exist with sources ju mping; between them 
and having different correlation slop es l|Russell et al.ll2007bl : 
ICoriat et al]|201ll ; [Gallo et alj|2012h . which implies that the 
jet luminosity is affected by additional, as yet undefined pa- 
rameters as well as the mass accretion rate and (possibly) 
the BH spin. The radio/X-ray (and IR/X-ray) correlation 
normalization therefore may not be an accurate method of 
estimating the relative ranking of compact jet power, how- 
ever the range in jet luminosities in the hard state (for a 
given X-ray luminosity) could still be improved. Neverthe- 
less, the results of FGR10 suggest there may not be a strong 
dependence of relative compact jet power on BH spin, even 
given the above caveats. 

It was also predicted by iBlandford fe Znaiekl (|l977l ) 
that jet power also scales with the square of the magnetic 
field strength; Pj ct oc a,B 2 , because the mechanism relies 
on differential rotation between a spinning BH and its mag- 
netosphere. One would therefore not expect a Pj ct oc a'i 
correlation to exist unless the magnetic field is nearly con- 
stant across all sourc es (this caveat is a lso mentioned, but 
then later ignored, in ISteiner et alll2013l ). In the hard state, 
it may be that B 2 oc rh, the mass accretion rate (e.g . 
iHeinz fc Sunvaevl 120031 ), and in ISikora fc Begelmanl 
it is argued that jet power is mostly driven by the magnetic 
flux, not the spin or accretion rate. Jet luminosity (measured 
in radio or IR) does scale positively with X-ray luminosity 
and hence with mass accretion rate, in the hard state. How- 
ever this is no indication of the Blandford & Znajek mech- 
anism, as a positive relation between X-ray luminosity and 
jet power is of course also expected for the Blandford-Payne 
mechanism (in which accretion powers the jet). In addition, 
iNixon fc Kind (|2013T ) have recently argued that jets that 
precess on short timescales cannot be spin powered. 

Even if BHXB jets were powered by BH spin, the com- 
bination of several effects - including but not limited to 
uncertainties in measuring the relevant parameters - could 
still prevent us from finding a positive correlation at present. 
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